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Abstract 


We  show  that  any  conservative  scheme  for  solving 
scalar  conservation  laws  in  two  space  dimensions  which 
is  total  variation  diminishing  is  at  most  first  order 
accurate . 


IV 


1 .   Introduction 

We  study  finite  difference  schemes  for  computinq  weak 
solutions  (e.g.  solutions  with  shocks)  to  scalar  conservation 
laws  of  the  form 

(1.1)         ut  +  f(u)x  +  g(u)   =  0 

in  two  space  dimensions.   In  one  dimension  Harten  [  3  ]  has  given 
conditions  which  guarantee  that  a  scheme  is  strongly  nonlinearly 
stable  in  the  sense  that  the  total  variation  of  the  numerical 
solution  does  not  increase  in  time.   He  has  also  shown  how  to 
construct  schemes  which  have  this  property  and  are  second  order 
accurate  on  smooth  solutions. 

By  contrast,  in  two  dimensions  such  schemes  do  not  exist. 
We  show  that  any  conservative  scheme  which  is  total  variation 
diminishing  (TVD)  is  at  most  first  order  accurate.   Briefly,  the 
proof  goes  as  follows.   For  any  two  dimensional  scheme  there  is 
an  associated  one  dimensional  scheme  with  the  same  order  of 
accuracy.   If  the  2D  scheme  is  TVD  then  the  corresponding  ID 
scheme  is  monotone  at  least  on  certain  initial  data  and  therefore 
at  most  first  order  accurate  by  the  results  of  Harten,  Hyman  and 
Lax  [4] . 

While  it  is  not  logically  necessary  for  a  scheme  to  be  TVD 
in  order  to  be  total  variation  stable,  or  to  converge,  we  know 
of  no  method  for  computing  weak  solutions  that  is  not  TVD  and 
yet  can  be  shown  to  converge.   In  any  case,  diminishing  total 
variation  is  a  very  desirable  property  which,  needless  to  say, 
is  shared  by  the  exact  solution  (see  the  appendix  by  Keyfitz 
in  [4])  . 


In  Section  2  we  review  and  slightly  extend  some  known 
results  on  monotone  schemes.   These  are  used  in  Section  3  to 
prove  our  main  result.   We  also  make  some  remarks  on  its 
significance  and  limitations. 

2 .   Monotone  Schemes 

A  finite  difference  scheme  for  the  problem  (1.1)  is  said 
to  be  in  conservation  form  if  it  is  of  the  form 


,  ^    -.  ^  .   n  +  1     n      ,  ,^-n        ^-n   .    ,   ,— n       — n   , 

{2-l)  Uj,k  =  Uj,k  "  Xx(fj+l,k  "  fj,k>  "  V^k  +  l'^k' 


where 

X   =  At/Ax,   X   =  At/Ay 
X  y         J 

and  the  numerical  fluxes  f  and  g  are  functions  of  neiqhborinq 
values  of  u: 


(2.2)     f •  i  =  f ( n •    i     t  • • •  i    v  ■  ,   i  , 
],k      L-]-p,k-q        u-]+r,k+! 


— n    _  — .  n  n 

gj,k  :  gl  j-p,k-q  '•••'  uj+r,k+s' 


We  assume  f  and  g  are  consistent  with  f  and  q, 

(2.3)     f(u,...,u)  =  f(u),  g(u,...,u)  =  g(u)  , 
and  are  Lipschitz  continuous  in  all  arguments. 


The    discrete    L      norm    is    defined    by 
1 


(2.4  a)  II U  II  =    Ax    Ay       £ 

Ll  j,k 


jk 


or    by 
(2.4b) 


III! 


Ax    I     |U 


in  one  space  dimension.   The  total  variation  of  the  qrid 


function  U ,  a  discrete  approximation  to 
I Vu I  =  I u    +  lu  I,  is  given  by 

II  I     j,  I  I     y  I  J 


7u I dx  dy  with 


(2.5)  TV(U)     =     J       [Ay|Uj+lfk    -    Ujfk|     +    Ax|Uj;k+1    -    U. ^ 


or  in  one  dimension  by 


TV(U)  =  I     |Uj+1-U. 


We  denote  the  scheme  (2.1)  symbolically  by  u     =  F(u  ) 
Since  we  are  studying  a  single  timestep  we  will  generally 
suppress  the  index  n  and  use  the  notation 

U  -  F(u)   • 


We  recall  two  definitions. 

Def  2.1.   The  scheme  (2.1)  is  monotone  if  each  component 
of  u  is  a  monotone  nondecreasing  function  of  \j,    i.e.,  if 


2-6)       uj,k±  vj,k     v^'k  -  uj,k  1  vj,k      v^'k  • 


Def  2.2.   The  scheme  (2.1)  is  L.  contracting  if  for  all  U 


and  V, 


(2.7)       IU-VI    <  IU-VIT 

Ll  _       Ll 


These  two  properties  are  equivalent. 

Prop.  2.1.   The  scheme  (2.1)  is  monotone  if  and  only  if  it 
is  L,  contracting. 

Remark.   This  proposition,  and  its  proof,  are  closely  related 

to  the  lemma  of  Crandall  and  Tartar  [2]  which  appears,  for  example, 

in  [1]  . 

Proof .   Assume  that  lU-Vll    <  <»  ,  for  otherwise  (2.7)  is 

Ll 
automatically  satisfied.   It  suffices  to  consider  the  case  when 

U  and  V  differ  at  only  a  single  point.   Monotonicity  on  general 

data  is  equivalent  to  monotonicity  on  such  data  and  similarly 

for  L,  contraction.   This  follows  by  induction  on  the  number  of 

points  at  which  U  and  V  differ.   One  must  use  the  fact  that  if 

W  differs  from  U  at  one  set  of  mesh  points  and  W  differs  from  V 

at  a  disjoint  set  of  meshpoints  then  B  U— VU  T   =  I U-WH    +  Bw-vH   . 

Ll        Ll        Ll 

We  suppose  U.,  =  V..  for  (j,k)  f    (0,0)  and  that 

(2-8)         uoo  >  voo  • 

From  the  conservation  form  (2.1)  of  the  difference  scheme  it  follows 
that 

,2-91      ^k(V  "  V  =  ^k  (V-V>  =uoo-voo  • 

Using  this  we  obtain 


iiu-v  ii 


=    y    iu .,  -  v .. 

l       j,k       ^  ] 


>  |  I      (u -k  -  v -k) 

(2.10)  "    j,k         =>*  ^ 


uoo  ~  voo 


i.e. 


(2.11)  llu-VilT       >     IIU-V  U  T 

Ll  _       Ll 


with  equality  holding  if  and  only  if  U  .   -  V    are  all  of  the 
same  sign  (which  must  then  be  positive  by  (2.8)  and  (2.9)). 
Comparing  (2.7)  and  (2.11)  we  conclude  that  the  scheme  is  L.. 
contracting  if  and  only  if  it  is  monotone,  completing  the  proof, 

In  Section  3  we  will  need  to  consider  schemes  in  one 
dimension  of  the  form 


(2.12)        Un+1  =  Un  -  \     (f?,,-  fn) 
3       3     *   3+1    3 


which  have  a  weaker  monotonicity  property. 

Def  2.3.   A  scheme  of  the  form  (2.12)  is  monotone  on  monotone 
data   (MMD)  if  it  is  monotone  in  the  sense  of  (2.6)  for 
all  monotone  data  U  and  V,  i.e.  if 


U  .  ,  >  U  .  ,   V  .  ,  >  V  .  ,   U  .  >  V  .     v. 

D+l  -   ]'    j+l-j'    ]-]      ] 


implies 


U  .  >  V  .    Vj 
D  -   3 


The  proof  of  Proposition  2.1  shows  that  if  (2.12)  is 
L,  contracting  on  monotone  data,  then  it  is  RMD.   Combinina 
this  with  the  results  of  Harten,  Hyman ,  and  Lax  [  4  ]  gives 

Prop.  2.2.   A  scheme  (2.12)  which  is  L,  contracting  on 

monotone  data  is  at  most  first  order  accurate  at  any  point 

where  U.    is  a  twice  dif f erentiable  function  of  U  . 
3 

Remark.   There  are  many  schemes,  such  as  the  original 
Godunov  scheme,  which  have  Lipschitz  continuous  fluxes  which 
are  not  continuously  dif f erentiable  everywhere.   However,  for 
all  known  schemes,  the  fluxes  are  infinitely  dif ferentiable 
except  at  a  few  points  (e.g.  sonic  points) .   Consequently  we 
believe  there  are  no  second  order  accurate  schemes  which  are 
L,  contracting  on  monotone  data. 

There  may  be  an  alternative  proof  of  this  which  does  not 
depend  on  differentiability.   We  conjecture  that  monotone 
schemes  exhibit  false  entropy  production  for  smooth  solutions 
which  is  0(h)  whenever  the  fluxes  are  Lipschitz  continuous, 
thus  implying  first  order  accuracy.   (see  Crandall,  Majda  [1] 
for  entropy  fluxes  in  2D) .   We  do  not  pursue  this  approach 
here  but  rather  follow  the  proof  of  Harten,  Hyman  and  Lax  [  4  ] , 
who  show  that  a  scheme  which  is  L,  contracting  on  all  data  is 
at  most  first  order  accurate  (again  at  points  of  differentiability, 
though  they  don't  make  an  issue  of  this). 


Proof:   If 


3U  . 
(2.13)       ^=p-   >    0      for  all  j  ,k 
k 

when   U,  =  Urt  is  constant  in  k  then  by  [  4  1  the  truncation 
k     0 

2 
error  is  of  the  form   h(B(u)u  )   +  0(h  )  where  B  ^  0  is  an 

explicit  continuous  combination  of  partials  of  the  form  (2.13) 

A  scheme  which  is  monotone  on  monotone  data  satisfies  (2.13) 

whenever  u  is  monotone.   This  gives  (2.13)  when  U  is  constant 

by  continuity  of  the  derivatives  since  constants  are  limits 

of  monotone  functions.   This  completes  the  proof. 


3 .   First  order  accuracy  of  TVD  schemes  in  2D. 

We  state  the  main  result  imprecisely  to  avoid  a  complicated 
(and  unenlightening)  statement. 

Assertion  3.1.   Except  in  certain  unlikely  cases, 
any  TVD  scheme  of  the  form  (2.1)  in  two  dimensions  is  at  most 
first  order  accurate. 

Proof.   From  any  scheme  of  the  form  (2.1)  for  the  2D 
problem  (1.1)  we  can  construct  a  scheme  of  the  form  (2.12)  for 
the  ID  problem   u   +  f (u)   =  0   that  has  the  same  order  of 
accuracy.   We  show  that  if  the  2D  scheme  is  TVD  then  the 
resulting  ID  scheme  will  be  L,  contracting   on  monotone  data. 
This,  together  with  Proposition  2.2,  proves  the  assertion. 


Suppose  we  apply  (2.1)  to  a  qrid  function  U  which  is 
independent  of  y, 


(3.  1)  n    =  V.      Vj  ,k 

jk     3 


Then  (2.1)  reduces  to  the  following  ID  scheme  for  V: 


(3.2)       V.  =  V.  -  X  (f  ...-f .) 
3     3     xv  ]+l  y 


where 


f  .  =  f (V.   , . . . ,V. .) 

3       3-P       3+r 


f  (U.    ,    ,  .  .  .  ,U.    .  _,)  . 
D~p,k-q'      j+r,k+s 


The  value  of  k  chosen  in  evaluating  f  is  arbitrary  by  (3.1). 
From  (2.3)  we  see  that  f  is  consistent  with  the  flux  function  f. 

f  (u,  .  .  .  ,u)  =  f  (u)  . 

Moreover  (3.2)  will  be  a  second  order  accurate  scheme  for 
u   +  f(u)   =  0  if   (2.1)  is    second  order  on  (1.1). 

The  total  variation  as  defined  in  (2.5)  is  the  sum  of  the 
x  and  y  variations, 

(3.3a)       TVx(u)  =  Ay  J       |Uj+lfk  -  Uj#k| 

3  » K 

and 

(3.3b)       TV  (U)  =  Ax   T   |U.  ...  -  U.  ,  |  =  7  IIU   ...  -  U   JT 
y  j,k    : '       l  D  '      k    '  '        '  '    1 


where  II  •  II    is  the  ID  L,  norm  (2.4b). 
Ll  l 

We  construct  data  for  which  TV   is  conserved 

x 


Then 


by  (3.3b)  total  variation  in  2D  is  related  to  L..  contraction 
in  ID.   Suppose  V  and  W  are  two  one-dimensional  grid  functions 


satisfying 


(3.4 


V.  =  W.  =  u 
J  J 


V.  <  V.    ,  W.  <  W..., 

:  -  D+i   :  -  ]+i 

v.  =  w.  =  un 

V.  >  V.  ,.  ,  W.  >  W.  ., 


V.  =  W.  =  u 

3     D 


for   j  <  -2N  , 


for   -2N  <  j  <  -N  , 


for   -N  <  j  <  N  , 


for   N  <    j  <    2N  , 


for   i  >  2N 


where   u0  >  u   and  N  is  large  relative  to  the  stencil  size 

U        oo  J 

determined  by  p,q,r,  and  s  in  (2.2)  .   Suppose  that  L  and  M 
are  also  large  integers,  and  define  U  from  V  and  W  to  be 
constant  in  y  along  each  of  2M  horizontal  strips  of  width  L, 
taking  the  values  of  V  and  W  alternately 


(3.5) 


V.   if   k  =  2k, L  +  k_ 
D  1      2 


with   0£k   <  M,  0  <_   k   <  L, 


W.   if   k  =  (2k1+l)L  +  k2 

with  0  £  k,  <  M,  0  <_  k2  <  L, 


u    otherwise 


Then 


II  U-u  II .  < 


and 


(3.6a)       TVx(li)  =  4LM (Ug-uJ Ay 


(3.6b)        TV,  (U)  =  (2M-1)  "V-W'l    +  "  V-u  J1  T  +  "w-u^1. 

y  lx        ^       l1 


Here  we  mean  the  one-dimensional  L,  norm  (2.4b). 

Since  N  is  large  compared  to  the  stencil  size  determined 

by  p,q,r,  and  s,  and  U.,  is  constant  in  j  and  k  for 

3K 

| j |  <    N  and  0  £  k  £  2LM  it  follows  from  the  consistency  of 
(2.1)  that 


U    =  u   for   q  <  k  <  2LM-S 


Since  U.,  =  u   for  I j I  sufficiently  large  it  follows  that  the 

jk     °°      '  J  '  J 

one-dimensional  variation  satisfies 


TV(U   ,  )  >  TV(U   .  )   for   q  <  k  <  2LM-S 
and  hence 


(3.7; 


TV  (U)  >  TV  (U)  -  C 
x     —    x 


]  0 


where  the  constant  C  depends  only  on  possible  edge  effects  near 
k  =  0  and  k  =  2LM  and  is  independent  of  L  and  M. 

In  the  y-direction  we  obtain  a  lower  bound  on  the  total 
variation  of  u  by  comparing  values  at  the  middle  of  successive 
strips  where  the  2D  scheme  reduces  to  the  ID  scheme  (3.2)  for 
L  sufficiently  large.   Suppose  V  and  W  result  from  applying 
the  ID  scheme  (3.2)  to  V  and  W  respectively.   Then  we  obtain 


(3.8)  TV  (U)  ^  (2M-1)  II  V-Wll 

y  l1 

Now  suppose  the  ID  scheme  is  not  L.  contractina  on  V  and  W  and 


that  consequently 


(3.9)  HV-WllT       =    II  V-Wll  _       +    e 

Ll         Ll 


with   e  >  0.    Then  by  (3.6b),  (3.7),  (3.8)  and  (3.9)  , 


TV(U)  =  TV  (U)  +  TV  (U 
x         y 


>  TV  (U)  +  TV  (U)  +  2(M-l)e  -  C, 
—    x         y  1 


where  C^  =  C  +  Hv-uJ1  +  II W  -  ujl  .    By  taking  M  sufficiently 
large  we  obtain 

(3.10)        TV(U)  >  TV(U) 

If  the  2D  scheme  is  TVD  then  (3.10)  cannot  occur  and  hence  the 

corresponding  ID  scheme  must  be  L,  contracting  on  all  data  V,W 

of  the  form  (3.4)  .   Since  we  can  take  V.  =  W.  for  j  <  0  or  for 


11 


j  >  0, in  particular  the  scheme  must  be  L,  contracting  on 
monotone  data.   By  Proposition  2.2  the  ID  scheme  (and  hence 
the  2D  scheme)  is  at  most  first  order  accurate. 

Remarks. 

(1)  A.  Harten  has  pointed  out  to  us  that  our  restriction 
to  explicit  schemes  is  probably  not  essential,  since  implicit 
schemes  are  also  rather  "local".   He  conjectures  that  our 
proof  would  apply  to  implicit  schemes  as  well  by  taking  L ,  the 
strip  width,  to  be  very  large. 

(2)  Despite  our  results,  split  schemes  in  2D   seem  to  work 
quite  well  for  gas  dynamics.   It  would  therefore  be  interesting 
to  prove  stability  and  convergence  for  some  second  order  accurate 
2D  scheme  for  scalar  conservation  laws.   This  was  in  fact  our 
original  intention. 

(3)  The  differences  between  scalars  and  systems  are  more 

pronounced  in  higher  dimensions  than  in  one  dimension.   For 

2 
example,  the  Euler  equations  are  not  well  posed  in  L,  (IR  ),  and 

gas  dynamical  shock  fronts  are  more  stable  than  scalar  shock 

fronts  in  two  dimensions  [5]. 
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